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Abstract 

Using a data sample with integrated luminosity of 3.9 fb _1 collected in e + e~ 
annihilation with the CLEO-II detector at the Cornell Electron Storage Ring, 
we have measured the branching ratios for the decay modes Df — > (rj, rj')ir + 
and Df — > (rj, r/')p + relative to Df — > 4>ir + . These decay modes are among the 
most common hadronic decays of the D+, and can be related by factorization 
to the semileptonic decays Dj — > (r],r]')£ + ui. The results obtained are com- 
pared with previous CLEO results and with the branching ratios measured 
for the related semileptonic decays. We also report results on the Cabibbo- 
suppressed decays of the D + to the same final states. 
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I. INTRODUCTION 



Among the most common hadronic decay modes for the are the decays — ► 
(7],r]')ii + and Df — > (r),r)')p + , where the notation — ► (r],r] r )ir + represents the decays 
-Dj" ~~ ^ r ? 7r+ and -D+ — > r/'7r + . As can be seen from Figure |l], they are related by factorization 
to the semileptonic decays Df — > (77, r/)£ + z//. This relation has been extensively discussed 
by Kamal, Xu, and Czarnecki 0. One prediction of the factorization hypothesis is that the 
D+ decay rate to r/p + can be simply related to the corresponding semileptonic decay rate 
evaluated at q 2 = m 2 : 

r p+ _> = 6^alf 2 p \V ud \ 2 - ^)l(^)- (1) 

Here / p is the decay constant for the p and a% is a strong interaction coefficient that is 
measured in two-body hadronic D° decays. 



3280397-005 



D 





FIG. 1. Related Feynman digrams for hadronic and semileptonic Df decays 

To test this factorization prediction experimentally, a shape for the form factor must be 
assumed. It is expected to be very close to the form factor for D° — > K~e + u e , for which 



T(D C 



K'e + u P 



dT(D° -> K-e+v e )/dq 2 



1.30 ±0.01 GeV . 
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This number is calculated using the CLEO measurement of the form factor ||. Assuming 
a similar pole shape for the form factor yields the prediction T(Df — > f]p + )/T(Df — > 
r]e + u e ) « 2.9 and r(L>+ -> rfp + )/T(D+ -> r/e+z/ e ) « 2.9 §. 

In 1992 CLEO Q measured the branching ratios for the hadronic modes studied here 
using a much smaller data sample of 0.69 fb _1 . Combining these measurements with the 
more recent CLEO measurements of the semileptonic modes ||, we calculate the T(Df — > 
T]p+)/V(D+ -> r/e+i/ e ) = 4.3±1.1 and T(D+ -> rfp + )/T(D+ -> »7'e + i/ e ) = 14.8±5.8. The last 
number is well above the factorization prediction of 2.9. Theoretical efforts to understand 
charm nonleptonic decays, even including final state interactions, cannot account for the 
very large branching ratio for — > r]'p + p,ID] . 

Because of the interest in these branching fractions, we have remeasured them using the 
much larger data sample now available. In the present analysis we use data corresponding 
to an integrated luminosity of 3.9 fb _1 (which includes the 0.69 fb~ used in the previous 
analysis) to remeasure the four modes, D+ — > (r],7]')n + and Df — > (r],T)')p + . The data 
were collected with the CLEO II detector at the Cornell Electron Storage Ring (CESR), at 
center-of-mass energies equal to the mass of the T(4S) and in the continuum just below the 
T(4S) resonance. 

The CLEO-II detector is designed to detect both charged and neutral particles with 
high resolution and efficiency The detector consists of a charged-particle tracking system 
surrounded by a time-of-flight scintillator system. These are in turn surrounded by an 
electromagnetic calorimeter which consists of 7800 thallium-doped Csl crystals. This inner 
detector is immersed in a 1.5 T solenoidal magnetic field generated by a superconducting 
coil. Muon detection is achieved using proportional tubes interleaved with iron. A more 
complete description of the detector can be found elsewhere H. 



II. EVENT SELECTION 

All events in this analysis are required to pass standard CLEO criteria for hadronic events. 
Since all the signal modes involve only pions in the final state, systematic errors are reduced 
by imposing no hadron identification cuts on either signal modes or the normalization mode. 
All candidates are required to have x = PD S /P m ax > 0.63 (P^ax = ^beam ~~ ^d s ) t° 
suppress combinatoric background. Throughout this paper, reference to a particular charge 
state implies the inclusion of the charge-conjugate state as well. 

All photons are required to be in the good barrel region of the calorimeter (| cos 9\ < 0.71), 
to have a minimum energy of 30 MeV, and to not match the projection of a charged track. 
We choose pairs of photons whose invariant mass is within 2.5 a(M) of the nominal tt° mass; 
a(M) is approximately 6 MeV/c 2 . We then kinematically constrain the 77 pairs to the 
nominal tt° mass in order to improve the momentum resolution of the 7r° . We also require 
that I cos^ol < 0.8, where 8 n o is the angle between one 7 in the 7r° rest frame and the ir° 
momentum in the lab frame. The signal is flat in | cos^o| and the background peaks toward 
+1. 

For i] — > 77 decays, the t] is selected in a manner similar to the ir°, but with the additional 
constraint that photons which could be paired to make 7r°'s with momentum greater than 
0.8 GeV/c are rejected. We also detect 77's using the 77 — > 7r + 7r~7r decay chain, although 
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this mode gives a sample with fewer events and less significance than the two-photon decay 
mode. A tt° momentum greater than 0.4 GeV/c is required. All 77 candidates within 2.5 
a(M) of the nominal mass are considered, where a(M) is the r.m.s. mass resolution for the 
given mode, typically about 14 MeV/c 2 for the 77 mode, and 6 MeV/c 2 for the 7r + 7r~7r° 
mode. In order to improve the momentum resolution of the 77, the decay particles from the 
77 are kinematically constrained to the nominal 77 mass. 

To select rj' candidates we use the r]iT + 7i~ final state, where the 77 is detected in both 77 
and 7r + 7r _ 7r modes. Both 77 and 77' candidates are kinematically constrained to the nominal 
mass in order to improve the momentum resolution. 

Reconstruction efficiencies and invariant mass resolutions were determined by using a 
GEANT-based [0 Monte Carlo (MC) simulation of the detector. 

III. D+ DECAYS INTO MODES CONTAINING A tt+ 

Five modes are studied in which a pion is produced in the weak decay: 



1. 


Dt- 


-* 07T + 


(the 


normalization 


mode), - 




2. 


Dt~ 


-> 7]7l + 


> v - 


* 77 






3. 


Df~ 


-> 777T + 


, v - 


> 7l + 7l~7l° 
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-> 1]'7l^ 


', 7]' - 


-> 777r + 7r _ , 77 — 


> 77 




5. 


Dt- 




', 77' - 


-> 777r + 7r _ , 77 — 


> 7T + 7r _ 7r 





We require the pions that come directly from the weak decay to have momentum greater 
than 0.7 GeV/c and the 77 or 77' from the to have momentum greater than 1 GeV/c. This 
reduces the background from random combinations. 

A. D+ -> 

Since this decay involves a pseudoscalar meson decaying into a vector meson and a 
pseudoscalar 7r + , the must be polarized in the helicity zero state. We take advantage 
of this by cutting on cos 6^+, where 8k+ is the angle between the K + momentum and the 
direction opposite to the momentum in the rest frame. The angle is shown in Figure |[ 
The signal has a cos 2 9 K + distribution, while the background is flat in cos 9 K + . We require 
I cos6 K +\ > 0.45. 

We select <fi mesons within ±8 MeV of the peak mass, and form the <f)ir + mass spectrum 
shown in Figure The 07r + mass distribution shows two clear peaks, one from the and 
the other from the D + . To fit the spectrum we use four functions: 

1. The -D+ signal is fit to a sum of two Gaussians with a common mean; the widths and 
relative areas are fixed to values determined from the Monte-Carlo signal simulation. 
The mean is allowed to vary in the fit. 
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FIG. 2. Illustration of the helicity angle, 9k+- All vectors represent momenta in the 4> rest 
frame. 

2. The D + signal shape is of the same form as for the Df, with the mass constrained to be 
0.099 GeV/c 2 less than the mass, which is the precisely measured mass difference 

3. The shape of the function used to represent the — > (<fi,r),r)')p + feed-through is 
determined from Monte Carlo. This feed-through causes a broad peak in the mass of 
the (4>, rj, rj')ii + system centered at 1.7 GeV/c 2 , which is parameterized with a Gaussian. 
The normalization of the feed-through is determined from the measurement of the 
branching ratio [|J. 

4. A second-order Chebyshev polynomial is used to represent the combinatoric back- 
ground. 

This fit yields 3748 ± 91 events. In all other fits, four functions are also used, although 
the combinatoric background shape depends on the particular mode. 

B. D+ — >■ TJ7T+ 

In Figure f| we show the r]TT + invariant mass spectrum for both decay modes of the rj. 
The signal peaks are evident for both the Df and D + . The peak at the contains 766 ±44 
events for the channel rj — ► 77, and 154 ± 22 events for the channel 77 — ► 7r + 7r~7T°. Multiple 
entries into the plot from a single event are allowed, and no effort is made to select among 
them. In Table | we list the yields for different channels and their efficiencies for decay. 
We also list the measurement for the ratio T(D+ — > 7] 11 tt + )/T(D^ — > <f)ir + ). In the table e is 
the efficiency and eB is efficiency multiplied by branching fraction of the secondary decays. 
The systematic errors for the efficiencies relative to the 07r + mode have several sources 
and differ slightly from mode to mode. For the r] 77 7r + mode the systematic error includes 
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FIG. 3. The M{(p'K + ) distribution. The larger peak is due to the decay — > (J)tt + ; the smaller 
peak at lower mass is due to D + -> 

uncertainties in the relative charged track (4%) and photon detection efficiencies (5%). We 
studied the Monte Carlo shape by letting the width of the two Gaussians vary in the fit and 
then calculated the shift in the central value, giving us an uncertainty of 3%. We also used 
different background shapes to determine the uncertainty due to the unknown background 
shape, and obtained an error of 4%. The total systematic error obtained by adding these 
uncorrelated errors in quadrature is 8%. For the ^3 7r 7r + mode the systematic error includes 
uncertainties in the photon detection efficiency (5%) and in the signal (5%) and background 
(8%) shapes. In addition there was a systematic error of (10%) due to the modeling of 
multiple entries. The total systematic error obtained by adding these uncorrelated errors in 
quadrature is 15%. 



TABLE I. 


Fit Results for D+ - 


-> (r?,V> + - BR 


is the branching 


ratio of the 4>(r], rj') decay mode 


that is used. 










Mode 


N 


<%) 


(B{%) 


ryr(<H 


(f)7T 


3748 ± 91 


19.1 ±0.2 


9.4 




?7 77 7r 


766 ± 44 


9.6 ±0.1 


3.7 


0.52 ±0.03 ±0.04 




154 ±22 


4.5 ±0.1 


1.1 


0.35 ±0.05 ±0.06 


r/'(r? 77 )7r 


479 ± 26 


6.7 ±0.1 


1.1 


1.09 ±0.06 ±0.07 




58 ±9 


1.9 ±0.1 


0.2 


0.73 ±0.11 ±0.12 



The measured ratio for T(Df — > T]s 7r 7T + ) /T(Df — ► <p7i + ) shown in Table | is approximately 
two standard deviations lower than the corresponding ratio for the ?? 77 mode, taking into 
account the systematic errors which are not in common to the two modes. The measurements 
of — ► i]'tt + , described in the next section, show a similar discrepancy, as do the D + — > i]tt + 
and 7]'n + , although those have less statistical significance. As a result, we searched in some 
detail for a systematic discrepancy in reconstructing the two rj decay modes. To calibrate 
the relative efficiency for these modes, and to check the reconstruction program, we studied 
events of the type D* + — > D°ir + with D° — > K rj. This has a very large and clean signal, 
and an rj momentum spectrum very similar to that for the decays. Using this process, 
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FIG. 4. The M(r]ir + ) distribution for (a) r\ — > 77, (b) r\ — > 7r + 7r 7r°. 



TABLE II. Fit Results for D+ -» (rj,rj')n + . 

Mode N e(%) eB{%) r/r(07r) 

(fm 1133 ±72 20.3 ±0.2 9^9 

T7 77 7r 225 ±38 9.6 ±0.2 3.7 0.53 ± 0.09 ± 0.05 

r] 3n TT 50 ± 20 4.6 ±0.1 1.1 0.40 ± 0.15 ± 0.07 

r/'(r7 77 )7r 114 ±18 6.8 ±0.1 1.1 0.90 ± 0.14 ± 0.07 

7]'(V3n)Tr 12 ±7 1.9 ±0.1 0.2 0.52 ±0.29 ±0.09 



we measure B(rj —>■ 77)/!? (77 — ► 7r + 7r~7r°) = 1.53 ± 0.16 ± 0.10, compared to the PDG 
value of 1.64 ± 0.04 ||. This confirms that the relative efficiency for the two decay modes 
of the rj is reproduced properly in the Monte Carlo simulation. Other checks using the 
data also reproduced the expected ratio of B(rj — > jj)/B(rj — > 7r + 7r _ 7r°), although with 
limited statistical power. Since we were unable to isolate any systematic effect, we attribute 
the difference between the two — > i]tt + measurements to an unlikely set of statistical 
fluctuations. 

The yields and relative branching ratios for all of the D + decays into the same final states 
are shown in Table |Tj. The efficiencies for the D + modes are generally very close to those 
for the corresponding decays. 



C. D+ -> t/tt+ 

For this mode, we can apply cuts on both the rj mass and the rj' mass, reducing the 
background substantially. Each mass provides a kinematic constraint, helping to improve 
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the resolution for the 7]'n + mass. As a result, these modes are significantly cleaner than 
Df — > r]7T + . We require the momentum of the rj' to be greater than 1.0 GeV/c. 

In Figure [| we show the t/7t + invariant mass spectrum for both 77 decay modes. The peak 
at the Df mass contains 479 ± 26 events for the channel r] — > 77, and 58 ± 9 events for the 
channel 77 — > 7r + 7r~7r°. The efficiencies and relative branching ratios are shown in Table Q. 
The systematic error on the branching ratio measurement due to the uncertainty in charged 
track efficiency is negligible for the case of 77 — > 77 because the final state has the same 
number of charged tracks as the normalizing mode. The main contributions to the systematic 
error are the uncertainties in the photon detection efficiency (5%) and in the shapes used to 
describe signal (3%) and background (3%). The total systematic error obtained by adding 
these uncorrelated errors in quadrature is 6%. For the channel rj — > 7r + 7r~7r°, the main 
contributions to the systematic error are the uncertainties in the efficiency for charged tracks 
(4%) and photons (5%) and in the shapes for signal (10%) and background (4%), and in 
handling of multiple entries (10%). The total systematic error obtained by adding these 
uncorrelated errors in quadrature is 16%. 
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FIG. 5. The M{ri'iT + ) distribution, using the decay mode rf — » rjTr + ir with (a) rj — ► 77, (b) 
77 — > 7T + 7r~7r . 

In D s + — > i]'ti + , 7]' — > 7]7T + 7T~, 7] — > 7r + 7r~7T°, we found that there are many events with 
multiple combinations of pions which satisfy our selection criteria. Most of them come from 
real rj' decays in which different rearrangements of the same four charged pions (two directly 
from the 77' and two from the 17), plus the 7r°, pass our 77 and 7/ cuts. We take only one 
candidate per event, choosing the candidate with the minimum value of a x 2 based on the 
7T°, 77, and 77' masses: \ 2 = ( SM %' -f ( SM v) _|_ ( SM go) _ resu lting measurements are shown 
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in Tables | and 0. 



IV. D+ DECAYS INTO MODES CONTAINING A p+ 

The analogous Df decay channels, where the 7r + has been replaced by a p + , can be studied 
using very similar cuts. Because of lower rates, lower efficiency, and a serious problem 
with multiple combinations within the same event, the rj — > 7r + 7r~7r° decay does not add 
significantly to the measurements of these modes, and is not used. A data sample with 
about 20% less integrated luminosity was used for the measurements of these modes. 



A. D+ r]p+ 

For the decay mode — > r]p + , we need to consider the possibility of nonresonant tj-k+tt 
feedthrough. For Figure || we require the helicity angle to be in the range | cos 0^+1 > 0.45, 
and the invariant mass of the 7r + 7r° to be within ±170 MeV/c 2 of the p + mass. A fit to 
the resulting rjir + TT mass spectrum is shown, yielding 589 ± 43 D£ — > i]p + candidates and 
8 ± 32 D + — > 7]p + candidates; thus there is no evidence of D + — > r]p + . We cannot directly 
extract a branching ratio for — > r]p + , however, until we account for possible nonresonant 
feedthrough. 




FIG. 6. The M(r]p + ) distribution, with r] — > 77. 



Although cuts on the helicity angle and on the p mass region can be used, the most 
reliable way to measure the resonant branching ratio is to fit the Dalitz plot. By doing this 
we make full use of the di-pion mass and the helicity angle to isolate the r]p + signal. We 
therefore make a Dalitz plot of all events with 1.94 < M(r]7T + iT ) < 1.99 GeV/c 2 , removing 
the cuts on helicity angle and on the 7r + 7r° mass. In Figure |7| we show four Dalitz plots: 
(a) the signal region in the data, defined as 1.94 < M(r]ir + ir°) < 1.99 GeV/c 2 ; (b) the 
M Dg data sidebands, which are the mass regions 1.75 < M(r/7r + 7r°) < 1.90 GeV/c 2 and 
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2.04 < M(rjn + n°) < 2.24 GeV/c 2 ; (c) the full Monte Carlo of the r]p+ signal; and (d) 
a simulation using a parameterized Monte Carlo of nonresonant ?77r + 7r events generated 
according to phase space. 

We do not expect t]tt resonant structures in this Dalitz plot because isospin forbids 
ss — > rjir. For all four Dalitz plots, we recalculate the values of M 2 {r]Ti + ) and M 2 (7r + 7r°) so 
that the Dalitz boundary corresponds exactly to that of the mass of the [Q, giving the 
sidebands the same boundary as the signal region. This causes negligible smearing of the p + 
resonance. 

The most obvious feature of the Dalitz plot is that the p + region stands out so clearly in 
the data, even though there is a significant non-.D+ background which contains very little p + . 
A binned Dalitz fit to the data distribution in the signal region was performed using the sum 
of the distributions in the other three plots in Figure |7|. The normalization of the non-D+ 
component is fixed using a fit to the r]p + mass distribution as in Figure |]but without helicity 
angle and p mass cuts. The number of the resonant and nonresonant events is varied 



in the fit, with no interference term allowed. The results of the fit are shown in Table |HI . 
The systematic error includes uncertainties in the efficiencies for charged tracks (4%) and 
photons (10%) and the shapes for the signal (4%) and the background (3%). 

In order to understand the systematic error due to possible interference between the 
resonant and nonresonant decays we also did a coherent Dalitz fit. The density of the events 
in the Dalitz plot is represented by the expression 

I = A\ + A\ + B x 2AiA 2 cos(5i - S 2 ), 

where A\ and b\ are the amplitude and phase of the Breit-Wigner resonance, A 2 and 5 2 are 
the amplitude and phase of the nonresonant decay, both of which are assumed to be constant, 
and B is an additional constant which is allowed to vary from zero to one. The case B = 
corresponds to no interference between the resonant and the nonresonant part; the case 
B = 1 corresponds to full interference, expected if the nonresonant case were indeed a single 
partial wave with constant phase. The true case could lie anywhere between these two limits. 
In the fit when the constant B is allowed to float it takes the value 0.24 ± 0.20, consistent 
with no interference. We therefore use the result from the incoherent fit to determine the 
branching ratio, and use the result from the coherent fit with B = 0.44 to find a conservative 
systematic error from this source. This corresponds to a 3.6% error. The total systematic 
error, obtained by adding this error in quadrature with the other systematic errors mentioned 
above, is estimated to be 12%. 



B. D+ r]'p+ 

The decay — > r]'p + was reconstructed using the decay mode rf — > 7/7r + 7r~, with 
7] — > 77. We require the momentum of the rf to be greater than 1.0 GeV/c and the invariant 
mass of the two pions to be within ±170 MeV of the p + mass. In Figure |8| we can see a clear 
peak of Mrf^+^o. The fit yields 181 ± 18 D+ — > r]'p + events and —4 ± 10 D + — > i]'p + events; 
thus there is no evidence for D + — > r]'p + . 

As for the case of the i]p + decay mode, we need to subtract any nonresonant feedthrough 
into the r)'p + final state. In this case, however, a Dalitz plot is not as useful in separating 
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FIG. 7. Dalitz plot of Df — > 7]p + , with rj — ► 77. The horizontal axis is M% + ; the vertical axis 
is M^ n+ . (a) Data signal region; (b) Data Md s sidebands; (c) MC signal; (d) MC simulation of 
nonresonant 7]ir + ir°, generated according to phase space. 
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FIG. 8. The M(r]'p + ) distribution, with r( — > r]Tr + -K~ , rj — » 77. 

signal from background, because the kinematic range for the di-pion mass does not extend 
beyond the region of the p. We do not expect 77V resonant structures in this Dalitz plot 
because isospin forbids ss — > 7/71". We therefore fit the angular distribution alone to extract 
the p component. As for the case of the Dalitz fit for i]p + , we use three components in the 
fit: (a) the resonant signal shape, a fourth-order polynomial determined from the Monte 
Carlo, which includes the distortion of the pure cos 2 9 n + shape due to detector acceptance; 
(b) a nonresonant shape, which is linear; and (c) a non-Df background shape, which is 
a first-order polynomial determined by fitting the sidebands. As in the Dalitz fit, we fix the 
background normalization from the mass fit, and vary the normalizations of the signal 
and nonresonant parts. 

Figure |9| shows the fit of the helicity angle distribution for the events in the mass 
peak. The results of the fit are shown in Table |HI[ The total systematic error of 11% includes 
uncertainties in the photon detection efficiency (10%) and in the signal (4%), background 
(3%), and nonresonant (2%) shapes. 



TABLE III. Fit Results for D+ -> {r],rj')p + . 



Mode 


N 


e(%) 


eB(%) 


r/r(<M) 


(j)TT + 


3000 ± 81 


19.1 ±0.2 


9.4 




VP + 


447 ± 31 


1.82 ±0.07 


0.47 


2.98 ±0.20 ±0.39 


n' P + 


137 ±14 


1.36 ±0.04 


0.15 


2.78 ±0.28 ±0.30 



The yields and upper limits on the branching ratios for the D + decays into final states 
with a p + are shown in Table |TV| . 
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FIG. 9. The helicity angle distribution for events in the Df mass peak for the decay channel 
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TABLE IV. Fit Results for D + -> (rj, r]')p A 



Type 


N 


£(%) 


eB(%) 


r/r(</.vr)(90%aL) 


4>7T + 


970 ± 65 


20.3 ±0.2 


9.9 




VP + 


8 ±32 


2.1 ±0.1 


0.55 


<i.n 


rj ' p + 


-4± 10 


1.7 ±0.1 


0.19 


<0.86 



V. CONCLUSIONS 

We have measured with improved statistics the branching ratios of the two-body hadronic 
decays of the Df: — > r]7i + , r]'n + , r]p + , and r]'p + . The results are consistent with the 
previous CLEO measurements |4j] and have improved errors. Using weighted averages of the 
two 7] modes, our results for — > (r],i]')7c + are: 

0.48 ±0.03 ±0.04 



and 



and 

T(Df - 07T 



r (d+ - 


-> r]n + ) 


r(Df - 


-> 07T + ) 


iw - 


-> r]'ir+) 


r(Df - 


-> 07T+) 


» are: 




r W - 


-> r)p + ) 


iW - 


-> 07T + ) 


iW - 





1.03 ±0.06 ±0.07. 



2.98 ±0.20 ±0.39 



2.78 ±0.28 ±0.30. 
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Using these measurements and the published CLEO semileptonic measurements 0, we 
can calculate the ratios which test factorization: T(Df — > r]p + )/T(Df — > r]e + i> e ) = 4.4 ± 1.2 
and r(D+ — > r]'p + )/T(Df — > r/e + z/ e ) = 12.0 ± 4.3. The branching ratio for the mode — > 
f]'p + is much larger than the value of 2.9 expected from factorization. Using the normalization 
B(D+ -> 07T+) = (3.6 ± 0.9)% §, we calculate £(L>+ -> r/p+) = (10.0 ± 1.5 ± 2.5)%, where 
the second error is due to the uncertainty in the — > 07r + branching fraction. This 
branching fraction is very large, considering the flavor wave function of the rf is only partly 
ss and that the rate is suppressed for such a P-wave decay very close to threshold. There is 
no obvious mechanism by which final state interactions could cause such a large enhancement 
of one of the dominant decay modes. 

Table |V| summarizes the measurements of branching ratios for all four decays and 
compares them with theoretical calculations. Models which are successful in predicting other 
charm hadronic modes reasonably well predict B(D+ — > rj'p + ) to be 1 — 3% P~PH. This 
failure leads theorists to consider contributions to the amplitude from decay diagrams other 
than that shown in Figure 1. For example, Ball et al. |12| argue that the high branching 
ratio for — > rj'p could be due to a cs annihilation into a W + and two gluons, in which 
the two gluons hadronize as an rj' '. 



TABLE V. Measurements and Predictions for Branching Ratios of D^t Decays. 



Mode 


r/r(^+) 


BSW [11 1 


HK§ 


BLP 


//7T 
VP 

rj'ir 
rfp 


0.48 ± 0.05 
2.98 ± 0.44 
1.03 ±0.09 
2.78 ±0.41 


1.04 
1.96 
0.61 
0.55 


0.58 ±0.15 
2.86 ±0.71 
1.55 ±0.42 

u - 4O -0.32 


0.30 
1.83 
1.32 
0.59 



Using the normalization B(D + — > 07r + ) = (6.1 ± 0.6) x 10~ 3 @, we also calculate the 
D + branching fractions to the same final states. Table [VT] summarizes the results. Since the 
D + decays involve two diagrams which interfere, the theoretical calculations vary widely, 
and are expected to be somewhat less reliable than for the case. 



TABLE VI. Measurements and Predictions for Branching Fractions of D + Decays. The exper- 
imental upper limits are at the 90% confidence level. 



Mode 


Branching Fraction (%) 


BSW |11] 


HKi 


BLP 


//7T 


0.30 ±0.06 


0.004 


0.68 ±0.21 


0.34 


VP 


< 0.68 


0.06 


< 0.48 


0.01 


rj'ir 


0.50 ±0.10 


0.16 


< 0.48 


0.73 


v'p 


< 0.52 


0.05 


< 0.07 


0.12 
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